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ABSTRACT 


Based  on  the  theory  of  continued  fractions,  a  technique 
is  developed  for  the  reduction  of  high  order  multivariable 
systems.  The  mathematical  basis  for  which  these  techniques 
work  is  elucidated,  and  its  superiority  of  the  mixed  form 
over  any  other  form  of  continued  fractions  is  established. 

The  general  solution  to  linear  regulator  problem  is  developed 
and  the  properties  which  this  solution  exhibit  are  elucidated. 
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I .  INTRODUCTION 


The  capability  of  obtaining  simplified  mathematical 
models  for  use  in  the  analysis  of  high  order  dynamic  systems 
has  traditionally  relied  on  the  experience  and  ingenuity  of 
the  analyst.  Usually,  these  efforts  have  been  achieved 
using  both  frequency  and  time  domain  techniques. 

In  dealing  with  the  problem  stated  above,  and  through 
the  use  of  theory  of  continued  fraction,  this  research  has 
found  a  series  of  properties  applicable  to  reduction  of 
multivariable  systems,  lower  order  observers  and  derivations 
of  lower  order  systems  for  the  linear  regulator  problem. 

The  relationships  developed,  show  applications  in  areas 
where  the  high  order  systems  are  impractical  or  undesirable 
to  use  due  to  their  complexity  or  difficulty  of  implementation. 

In  Chapter  II  reduced  order  models  of  multivariable 
systems  for  the  first,  second  and  mixed  Cauer  forms  are 
developed.  Techniques  for  approximating  a  high  order  linear 
time-invariant  system  with  various  inputs  and  various  outputs 
by  a  reduced  order  model,  have  been  suggested  by  Chen  Cl], 
Meier,  L  and  Luenberger  [2],  L.  S.  Shieh  and  Y.  J  Wei  [3], 

M.  R.  Calfe  and  M.  Healey  [4],  L.  S.  Shieh,  J.  M.  Navarro 
and  R.  Yates  [5],  D.  A.  Wilson  [6],  L.  S.  Shieh  and  F.  F. 
Gaudiano  [7], 


9 


_  .j 


Most  of  these  methods  for  reducing  high  order  linear 
systems  are  based  on  the  following  principles: 

1.  The  low  performance  terms  can  be  discarded  and  the 
high  performance  terms  should  be  retained. 

2.  Linear  transformation  to  obtain  matrix  diagonaliza- 
tion  where  certain  diagonal  elements  can  be  dis¬ 
carded  . 

3.  The  sum  of  squares  of  the  errors  between  the 
responses  of  the  real  system  and  those  of  the 
approximate  model  at  the  sampling  instant  is 
minimized  in  order  to  obtain  the  parameters  of  the 
approximate  model. 

In  a  recent  paper,  Chen  [1]  proposed  a  reduction  of 
multivariable  control  systems  by  means  of  matrix  second 
Cauer  form  of  continued  fractions.  Through  the  method,  a 
simplified  model  is  obtained  by  keeping  the  first  several 
significant  matrix  quotients  and  discarding  the  others. 
However,  the  technique  (due  to  the  nature  of  the  Cauer 
second  form) ,  provides  satisfactory  results  in  the  steady 
state  region  only.  Furthermore,  M.  R.  Calfe  and  M.  Healey 
[4],  have  shown  that  the  method  does,  not  guarantee  the 
reduced  model  to  be  stable. 

In  Chapter  III  derivation  of  lower  order  system  for  the 
linear  regulator  problem  via  Cauer  form  is  obtained  and 
also  a  near  optimal  solution  for  the' original  system  can  be 
found  through  a  reduced  system. 
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II.  MULTIVARIABLE  SYSTEMS  REDUCTION  VIA  THE  CAUER  FORM 

A.  THREE  MATRIX  CAUER  FORMS 

L.S.  Shieh  and  F.F.  Gaudiano  [7]  have  shown  that  in  terms 
of  multivariable  systems  the  quotients  in  the  three  Cauer 
forms  are  replaced  by  matrix  quotients  and  the  division  in 
the  continued  fraction  is  replaced  by  matrix  inversion.  The 
first  matrix  Cauer  form  is 

T  ( s )  =  [H^s  +  [H’2  +  [H'3s  +  [H'4  +  C...]"1]"1]"1]"1]"1  (1) 

the  second  matrix  Cauer  form  is 

T(s)  =  [Hx  +  [H2|  +  CH3  +  [H4|  +  [. ..3"1]-1]"1]"1]"1  (2) 

and  the  mixed  matrix  Cauer  form  is 
1(3)  =  [Kx  ♦  qs  *  CKji  *  X'2  *  :k3  *  k'3s  *  k4| 

+  K»4  +  C...]"1]"1]'1]-1]"1  (3) 

where  ,  and  K’^  are  constant  m  x  m  matrix  quotients 

obtained  respectively  from  the  matrix  Routh’s  array  and  the 
generalized  Routh's  algorithm  shown  in  equations  (4a)  and 
(4b)  . 
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where  the  elements  of  the  first  and  second  rows  of  equations 
(4a)  and  (4b)  are  the  matrix  coefficients  of  the  system 
given  in  equation  (5). 

T(s)  =  CA2,n  Sn_1  4  A2n-lSn‘2  +...A23S2  +  A22S  +  A^lx 

(5) 

[A^,n+lSn  +  A1,nSn'1  +...A13S2  +  A12S  +  A11]"1 

The  elements  of  the  third,  fourth  and  subsequent  rows  in  (4a) 
and  (4b)  are  evaluated  respectively  for  the  three  Cauer  forms 
by  the  formulation  shown  in  (6a)  through  (6c). 

Hp  =  Ap,n+2-pAp'  +1,1  p  =  l,2,3,...n 

Aj,i  =  Aj-2,  i+1  -  Hj-2  Aj-l,i  j  =  3,4, ...2n 

det  Ap+1,1  ^  0  i  =  n+1,  n,  n-l,...l 

(6a) 


(6b) 


Hp  =  Ap ,  lAp"  +1,1 
Aj ,i  =  Aj-2 ,  i+1 
det  Ap+1,1  t  0 


p=l,2 ,3  ,  . . . 2n 
Hj -2Aj -1 ,  i+1  j  =  3  ,4 , .  . . 2n 
i-l,2,3,...n 


Kp  =  Ap,l  Ap+1,1 

K?p  =  Ap,  n+2  -  p  Ap+1,  n+l-p 

Aj  , i  =  Aj-2,  i+1  -  Kj -2A j -1 ,  i+l-KT  j-2  Aj  =  l,i 


p  =  l , 2 , 3  ,  .  .  .n 
j=3 ,4, . . .n+1 
i=l  ,2,3,  ... 


(6c) 


It  is  important  to  note  that  since  the  Cauer  first  and  second 
forms  are  special  cases  of  the  Cauer  mixed  form,  their 
formulation  in  (6a)  and  (6b)  can  be  derived  directly  from 
(6c),  by  letting  all  kp '  s  or  all  K£> '  s  equal  to  zero  respec¬ 
tively  . 

B.  STATE  SPACE  FORMULATION  FOR  THREE  CAUER  FORMS 

The  Cauer  Mixed  Form  -  Consider  a  typical  fe'edback  system 
with  a  minor  feed  forward  loop  as  shown  in  Figure  1.  The 
overall  transfer  function  is  given  by: 

T(s)  =  |||y  =  [G  +  F][I  +  (G  +  F)H]-1 .  (7) 

Equation  (7)  can  be  rewritten  as  a  mixed  matrix  of  continued 
fractions 

T(s)  =  [H  +  [F  +  G]'1]"1,  (8) 


14 


where 


t 


Figure  1.  Feedback  and  Feedforward  controls. 

If  the  subsystem  G  is  expanded  again,  equation  (3)  is  obtained. 
This  equation  can  be  represented  by  the  block  diagram  shown 
in  Figure  2.  Where  a  2  inputs-2  outputs  nth  order  system 
is  shown.  Again,  it  is  important  to  note  that  if  all  Ki ' s 
go  to  zero  in  Figure  2  the  block  diagram  representation  of 
Cauer  matrix  first  form  nth/ 2  order  system  as  shown  in 
Figure  3  will  automatically  be  obtained.  In  a  similar 
fashion  if  the  k'i's  go  to  zero,  the  block  diagram  represen¬ 
tation  of  Cauer  matrix  second  form  nth/ 2  order  system  as 
shown  in  Figure  4  is  obtained. 


Figure  2.  General  Matrix  representation  of  a  nth  order  system  by 
Cauer  mixed  form  with  two  inputs  two  outputs,  (n  even) 


Matrix 

second 


Going  back  to  Figure  2  and  allowing  e^  to  be  the  state 
variable  vector  (same  order  as  the  matrix  K  or  K*  ) ,  the 
equations  in  the  time  domain  can  be  written  as  follows: 


KlK2gl  +  KiK4"2  +  K1K6*3  +  Wn/2 


-(I  +  K,1K2  +  K1K’2)e1-(K,1K4  +  K.^ )  -  (K^Kg+K^g  )e3  -  ... 


-(K'lVKlKn)^n/2  -  KlK2el  "  KlK4e2- ‘  * -KlKnen/ 2  +  IU 


K’K*  e  +  (K1+K'  )K’  e„  +  OC+K'  )K’  e,+  .  .  .  +  (KI  +K1  )K'  e 
121  1342  1363  13nn/2 


-(K1K,2+K'1K2)e1-[(I  +  K,1+K,3)K4+(K1+k3)K'4]e2- 


[(Ki*K’3)K6H:<1«3)K'6]i3-....C(Ki.K'3)KnMK1tK3)K^en/2 


-K,  K,e.-(X,  ♦K->K„«,-(K,  *K,)K,-.  .  .-(K,  *K,)K  e  ,,*IU 
121  1342  136  13nn/2 


K^KV, e i+  ( + K  ^  )  K'4 e  2  +  ( K 'i + ^  +  IC5  )  K'g  e  3  +  .  .  .  +  ( K'i+ K’3 + K’5  )  K^en /  2 


-(K1K’2+K,1k-2)41-[(K'1+K,3)K4  +  (K1+K3)K[+Je2-[I  +  (K'1+K,3+K'5)K6  + 


(Kl+K3  +  K5)K6]43-*’--[(Kl+K3  +  K5)Kn+(Kl+K3  +  K5)K,n]®n/2 


(10) 


(11) 


'KlK2el_(Kl  +  K3)K4e2-(Kl+K3  +  K5)K6e3~ ’ ’ •_(K1+  K3  +  K5 5 Knen/ 2  +  IU 


(12) 
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KiK2®l+(Kl+K3 )^4®2+ (Kl+K3+K5 )K6®3+  * • • 


+  (K,l+K3+K5+-'-+Kk-l)Kken/2 


-(K'1K2+K1K'2)e1  -  C(K,1+K'3)K4+(K1+K3)K,4]e2 


[(K1+K3+K5)K6+(K1+K3+K5)K6]43-- • 


Cl  +  (K’+K'+K,,  +  .  .  ,+K*  ,  )K  +(K,  +K,+K,+.  .  .+K  .  )K’  ]e  /0 

loo  n-1  n  135  n-1  n  n/2 


-K1K2e1-(K1+K3)K4e2-(K1+K3+K5)Kge3- 


(K,  +K-  +  K-+  ...  +  K  JK  e  +  IU 
1  3  o  n-1  n  n/2 


Y  =  K2e1+K,2e1+K4e2+K'4e2+  . 


.  +K  e  /0  +  K*  e  /0 
n  n/2  n  n/2 


where  U  and  Y  are  the  input  and  output  variables  of  the 
system  respectively  and  I  is  the  Identity  Matrix.  From 
Equations  (10),  (11),  (12),  (13),  and  (14)  the  following 
corresponding  matrix  formulation  can  be  obtained: 

[A1]E  =  -[A2]E  -  [A3]E  +  [ 3 JU , 


Y  =  cc1j1e  +  [c2rE. 


where 


] 


KiK*2 

K^K’2  (KJ+K^)K^ 
K1K2  (K1+K3)K4 


KiK6 . K1KA 

<Ki+K3)K6 . (Kl+K3)Kn 

(K^K'+K«)K'...(K'1+K'3+K'5)K^ 


*1*2  (K'i+K')K-4(K^K*3  +  K-5)K'6.  .  .  (K^K'3+.  .  .-K^K^ 


(17) 


(18) 


I  +  K^K2  +  K1K,2 

K1K4  +  KiK4--* 

K1K2+K1K2 

i  +  (kl  +  k3)x’4  + 

(X*_  +  x'3)x 

K1K2+K1K2 

CKX  +  K3)K'4  +  ( K'i 

+  k'3)k4.  .  . 

K1K'2  +  K'1K2 

(Kx  +  K3)K’4  +  (K^ 

+  K'3)X  .  .  . 

•KlK6+KiK6*-- 

...(Ki+k3)k*+(k^k-3)k6... 

.  .  .I  +  (K1+K3  +  K5)K'6  +  (K’1+K'3  +  K'5)K6.  .  . 


.  (i<l+K3+K5  )k’6  +  (K'^+K3  +  K3  )Kg  .  .  . 
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.  .  .K,K'  +KLK 
In  In 

.  .  .  (KI  +K*  )K  +  (K,  +K,)K' 
1  3  n  1  3  n 


.  . . (Kt +K1+KV )K  + (K, +K,  +  KC 
1  3  5  n  13  5 


)  K 


T 

n 


.I+OC'+Kl  +K’+.  .  .+K'  ,  )K  +(K,+K,+K,+  .  .  .+K 
lob  n-1  n  1  3  o 


j 


(19) 


(20) 


$n/2 


1_ 

_j 

K1K2 

K1K4 

K1K6 

.  .K.K 

1  n 

K1K2 

(K1+K3)K4 

<WK6  •' 

. . (K,+K-)K 
l  3  n 

K1K2 

(Kx+K3)K4 

'WW- 

.  . (K, +K.+K,  )K 

1  3  5  n 

K.  Kr,  (K,+K-)K„  (K,  +K0+KC)KC _ (K.,+K„+Kc+. 


nr 

xi 

xm+l 

x2m+l 

x(n/2-l)m+l 

• 

X2 

* 

• 

42  = 

Xm+2 

*3  = 

x2m+  2 

en/2  = 

X(n/2-l)m+2 

X 

-  3L 

_X2m_ 

_x3m 

_x(n/2)m 

(27) 

*r 

Xm+1 

X2m+1 

x(n/2-l)m+l 

*2 

8  2  = 

Xm+  2 

e3  = 

X2m+2 

en/2  = 

x(n/2-l)m+2 

X 

_  m_ 

_X2m_ 

X, 

_  3m 

_x(n/2  )m 

(28) 

where  m  is  the  number  of  inputs  and  outputs.  It  is  important 
to  emphasize  the  properties  that  the  Cauer  matrix  mixed  form 
exhibits.  If  again  the  Ki's  go  to  zero  in  the  state  formula¬ 
tion  described  by  equation  (15)  and  (16),  they  will  reduce  to: 

[A  ]E  =  -[A213E  +  [B]U,  (29) 

and 

Y  =  [C21]TE.  (30) 

and  letting  K'i's  =  H'i's  equation  (29)  and  (30)  define  the 
state  space  formulation  of  the  Cauer  matrix  first  form  for 
an  nth/2  order  system,  where  n  is  even,  Equations  (29)  and 
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(30)  are  second  order  differential  equations  which  can  be 
simplified  to  first  order  differential  equations  by 
assigning  a  new  state  variable  [Z]  =  [El.  Thus  equations 
(29)  and  (30)  can  be  rewritten  as: 


[AX1]Z  =  CA21]Z  +  [B]U , 


Y  =  [C2i]  Z- 


where  [A^]  =  lim  [A^l  and  K'i's 
Ki’s  -*■  0 


=  H’i's 


then : 


H1H2  H1H,4 


HiH6 


H1H2  (Hi+H'3)H4  (Hi+H3)H6 
CAul-  H1H2  (H1+H')H'  (H^+H^+H^H’g. 


•  •  H',  H’ 

1  n 

.  .  (H‘  +H'  )H' 

1  3  n 

.  .  (H1,  +H'  +H,-)tf 
l  3  d  n 


H1H2  (H'^+H^  )H'l+  (H^+Hj+Hg  )Hg  , 


(H’,+H’  +HV+.  .  .H*  , 

135  n-1  n 


Z=  e. 


where  e^,  e^,  e^  •••  en/2  are  §iven  in  equation  (27). 

CA21]  =  lim  [A2]  =  til,  (35) 

Ki '  s-*-0 

where  I  is  the  Identity  matrix. 


(36) 

where  e^ ,  e^ ,  ...  e^2  are  given  in  equation  (28) 

[C01]T  =  lim[C0]T  and  K’i's  =  H'i's  [C0,  ]T  =  [H'  HV  H’H']. 
21  2  21  2  4  6  n 

Ki's-O  (37) 

This  formulation  given  by  equations  (30)  through  (37)  could 
have  been  obtained  directly  by  inspection  from  Figure  3. 

Similarly,  if  all  K’i's  go  to  zero  in  equations  (15)  and 
(16)  then 

E  =  -[A32]E  +  [ B ]U ,  (38) 

and 

Y  -  [C12]TE,  (39) 
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where 


=  lirnCA^]  and  Ki's  =  Hi's 
K'i's-M) 


H1H2 

H1H4 

. . .H. H 

1  n 

H1H2 

(h1+h2)h4 

(H].+H3  >H$ 

. . . (H, +H,)H 

1  3  n 

[A32]  = 

H1H2 

(H1+H3)H4 

<WH5 

)HC . . . (H1+H,+H,)H 
b  loon 

H1H2 

<H1+H3)H4 

“WH5 

)HC...(H,+H,+HC+...H  ,  )H 

6  13  5  n-1  n 

[Ci2]T  =  lim[C-^]^  and  Ki's  =  Hi's 
K'i's-MJ 

CC1  9]T  :  [H,  H.,  H,  .  .  .  H  ]  .  (41) 

1/  2  4  b  n 

which  define  the  state  space  formulation  for  the  Cauer  matrix 
second  form.  Again,  this  preceding  matrix  formulation  could 
have  been  derived  directly  by  inspection  from  the  block  dia¬ 
gram  shown  in  Figure  4.  Equations  (40)  and  (41)  have  the 
same  form  of  the  state  space  formulation  given  by  Chen  [1]. 

At  this  point,  the  following  observations  may  be  made  re¬ 
garding  the  state  space  formulation  just  developed: 

1.  The  elements  in  the  state  matrix  [Ai's]  are  simple 
matrix  combinations  of  the  matrix  quotients 
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obtained  from  the  continued  fraction  expansion  or 
Routh ' s  algorithm . 

2.  The  submatrices  appearing  below  the  main  diagonal 
have  the  same  value  as  the  submatrices  at  the 
diagonal  (with  exception  of  [Aj]). 

3  .  The  submatrices  which  appear  above  the  main  diagonal 
can  be  expressed  in  terms  of  matrix  quotients  in  a 
very  regular  way. 

4.  State  space  formulations  for  the  first  and  second 
Cauer  form  can  be  obtained  directly  from  the  Cauer 
mixed  form  through  direct  substitution. 

C.  APPROXIMATION  BY  THREE  CAUER  FORMS 

The  reduction  of  the  order  of  a  transfer  function  or 
decreasing  the  dimension  of  a  state  matrix  is  highly 
desirable  or  sometimes  necessary  in  the  analysis  and  design 
of  control  systems. 

In  terms  of  continued  fractions,  the  simplification 
problem  is  carried  out  by  expanding  a  given  transfer  function 
into  one  of  the  Three  Matrix  Cauer  Forms  of  continued  frac¬ 
tion  and  ignoring  some  matrix  quotients.  If  the  given 
system  is  outlined  in  state  variable  form,  the  simplifica¬ 
tion  method  is  realized  by  partitioning  the  matrix  and 
discarding  some  parts.  Several  examples  are  included 
for  demonstrating  the  power  of  the  method.  Also  a  thorough 
comparison  among  the  three  Cauer  forms  is  presented  and  their 
advantages  and  disadvantages  are  discussed. 
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Feedback  Gain  and  Feedforward  Gain  -  Consider  the  system 
shown  in  Figure  5.  The  closed  loop  of  the  overall  transfer 


function  is  known  as  follows: 
T(s)  =  [H1  +  [F2  +  G]"1]"1 


(42) 


Figure  5.  Block  diagram  for  a  typical  feedback  system 
with  a  minor  feedforward  loop  with  two  in¬ 
puts  and  two  outputs. 

where  G  =  0.  Equation  (42)  can  be  considered  as  the  simplest 
continued  matrix  fraction  expansion.  The  physical  meaning 
implied  in  the  formula  is  significant.  It  is  easily  seen 
that  when  F2  is  high  the  overall  gain  can  be  approximated  by 
H”1,  in  other  words,  dominates  the  behavior  of  che 

system. 
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This  fact  establishes  the  feedback  loop  as  the  most  im¬ 
portant  link  for  influencing  the  behavior  of  the  system, 
leaving  the  feedforward  loop  as  the  second  most  important 
link.  Furthermore  if  the  subsystem  G,  instead  of  being  zero 
is  still  a  high  order  transfer  function,  it  is  possible  to 
continue  the  expansion  one  after  another.  This  corresponds 
to  a  combination  of  many  feedback  and  feedforward  blocks  as 
shown  in  Figure  6.  It  should  be  noted  that  the  most  dominant 
term  is  and  the  second  influence  term  is  F2  .  When  the 
matrix  quotients  in  the  continued  fraction  are  lower  and 
lower  in  positions,  they  are  less  and  less  important  as  far 
as  the  influence  to  the  performance  of  the  system  is  con¬ 
cerned.  This  observation  is  the  general  basis  for  the  sim¬ 
plification  technique  developed  for  multivariable  systems. 

Considering  a  simple  case,  a  second  order  transfer 
function  such  as: 

T(s)  =  [A22S  +  A21][A13S2  +  Al2S  +  A11]~1  (43) 

can  be  expanded  into  three  different  matrix  Cauer  forms  of 
continued  fraction  as  follows: 


Cl(s) 

=  [H^S  • 

►  [h'2  + 

[H’3S  +  [H^]-1]-1]"1 

(44) 

C2<s) 

=  CHX  ♦ 

CH2i  ♦ 

CH3  +  [H^r1]-1]'1]"1 

(45) 

c3(s) 

=  CKX  ♦ 

K^S  + 

tK2i  *  kit1]-1 

(46) 
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where  F  ^(s),  Fg2(s),  and  Fc2(s)  are  the  first,  second  and 
mixed  Cauer  continued  matrix  fraction  expansions  of  equation 
(43)  . 


From  previous  observation,  it  is  known  that  the  most 
dominant  term  in  equations  (44),  (45),  and  (46)  is  H^s , 
and  K1  +  K'^s  respectively.  It  is  desirable  to  find  a 
meaningful  interpretation  for  the  dominant  term  of  each  one 
of  the  expansions  just  performed.  This  task  is  accomplished 
by  applying  the  initial  and  final  value  theorems  to  equations 
(44),  (45),  and  (46)  as  follows. 

1 .  Cauer  First  Form 

Performing  the  inverse  procedure  on  F  ,(s)  it  can  be 
written  as 


Fci ( s  5  =0Xsy=  CH2H3H4S+H2+H4  ]  Ch'iH2H3H4s2  +  (HiH2+H'iH4+H3H4 }  S+I 


(47) 


Applying  the  final  value  theorem  to  equation  (47)  and 
allowing 


U(s)  = 


1 

m 


where  m  is  the  number  of  inputs  is 
found  that 
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y  (t-*-°°)  =  [H  9+Hu  ] 
s-*-0 


(48) 


similarly,  applying  the  initial  value  theorem  and  allowing 


U(s) 


y(t-»0)  =  [H'  ]-1 
3^-00  L 


(49) 


The  meaning  implied  in  equations  (48)  and  (49)  is  very 
significant.  The  initial  conditions  dominate  the  behavior 
of  the  system.  In  other  words,  the  Cauer  first  form 
influences  very  heavily  the  transient  part  of  the  response. 

2 .  Cauer  Second  Form 

Performing  the  inverse  procedure  of  F^Cs),  this  can 
be  written  as: 
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Fc2 ( S 5 =UTiT= [H2+H4S+H2H3H43 [is2+ (H1H2+H1H4+H3H4 } S+H1H2H3H4]"1 

(50) 

Applying  the  final  value  theorem  to  equation  (50)  and  allowing 


n~\ 

n  n 

i 

1 

12 

—  it  is  found  that  y(t-*-«>)  =  [H. 
s  1 

*2 

. 

s->-0 

• 

1 

1 

m 

(51) 

Similarly,  applying  the  initial  value  theorem  and  allowing 


_11 

h 

h 

F2 

y(t^0)=CH2  +  H  ] 

S+00  ^ 

* 

• 

1 

1 

L"l 

__  m 

(52) 


The  results  obtained  in  equations  (51)  and  (52)  imply  that 
the  final  or  steady  state  value  dominates  the  behavior  of 
the  system.  In  other  words,  the  Cauer  second  form  influences 
very  heavily  the  steady  state  part  of  the  system  response. 

3 .  Cauer  Mixed  Form 

Performing  the  inverse  procedure  on  F  ,(s)  this  can 
be  written  as: 

Fc 3  ( S  )  =ijrfr  [K2S  +  K2  ]  [ ^1^2 S  2+  ( I  +  K1K2  +  K1K2  )  S+KlK 2  ]  ”1 


(53) 


Applying  the  final  and  initial  value  theorems  to  equation 
(53)  in  the  same  fashion  as  for  the  previous  forms,  the 
following  results  are  obtained: 


y(t-*-°o) 

s-*-0 


and 


y(  t-*-0 ) 

S  ->-oo 


(54) 


(55) 


These  results  show  the  steady  state  value  and  the  initial 
conditions  to  dominate  on  equal  levels  of  significance  the 
behavior  of  the  system.  Thus  the  Cauer  mixed  form  influences 
the  system  response  in  the  transient  part  as  well  as  in  the 
steady  state  part.  This  fact  makes  the  Cauer  mixed  form  a 
better  and  more  accurate  device  to  be  used  in  the  simplifi¬ 
cation  and  reduction  techniques  of  transfer  functions  of 
multivariable  systems. 


D.  TRUNCATION  OF  THE  CAUER  MIXED  FORM  THE  BEST  APPROXIMATION 
IN  A  DEFINITE  MATHEMATICAL  SENSE  FOR  MULTIVARIABLE 
SYSTEM  REDUCTION 

In  a  recent  paper,  M.J.  Goldman  and  C.T.  Leondes  [8] 
developed  the  mathematical  basis  for  the  simplification 
technique  involving  the  truncation  of  the  mixed  Cauer  form 
for  the  single  input  single  output  case.  In  this  section, 
an  extension  of  their  work  to  the  multivariable  case  is 
presented  and  its  superiority  over  any  other  form  of  con¬ 
tinued  fractions  is  established. 

The  transfer  function  matrix  T(s)  for  multivariable 
system  can  be  expressed  as  follows: 


T(s )  =  [A.  S11*1  +  A0  „  1Sn_2+.  .  .+A„,S2  +  A00S  +  A01  ]  x 

2  ,n  2 ,n-l  23  22  21 


+---A13S2  *  A12S  ♦  A^'1 


(56) 


where  Ai,j  are  constant,  m  by  m  matrices  and  Aij  =  aj[lj, 
j  =  l,2,...n+l,  where  each  aj  is  a  coefficient  of  the 
common-denominator  polynomial  or 

n+1 

A(s)  =  I  ajS2_:2  and  [I]  is  the  identity  matrix. 

3=1 

The  nth  convergent  of  a  mixed  matrix  Cauer  form  can  be 
represented  by  the  following  two  configurations: 
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(57) 


A  B^CX-.+X',  S  +  S[X„+X' S+S[ . . .S[X  +K' S]"1]"1]"1]"1 
n  n  11  i  i  nn 

A'  B'“1=i[K  i+K'  +i[X9i+K'0+i[  .  .  .i[K  -  +  X1  l-1]-1]-1]'1  (58) 

nnsl3ls2s2s  sns  n 


where 


A  =  X  A  1+SX’A  n+SA  ,,  ,  A  =  0  A  =  I 
”  n  n-1  ^  1  ~  1 


n 


n  n-1 


n-2 


B  =  X  B  .  +  SX'  B  ,  +  SB  0,  B  =0  B,  =  X,  +  X',  S 
n  n  n-1  n  n-1  n-2’  o  1  1  1 


(59) 


A'  = 
n 


— X  A'  ,  + 
s  n  n-1 


X’  A'  ,  +  —A  ,  A' 
n  n-1  s  n-2  o 


=  0  A',  = 


h 

s 


iK  B'  -1  *  K°n- 


n  s  n  n 


1  +  iBn-2’ 


B*  =  I  B' 


_  1 


— X  +  X'. 
si  ± 


(60) 


Recurrence  relations  (59)  and  (60)  have  been  derived  from 

standard  results  in  the  theory  of  continued  fractions 

(Rice  1964)  and  X  's  and  X'  Ts  from  the  generalized  matrix 

n  n 

Routh  Algorithm  (Shieh  and  Gaudiano  1974). 

Since 


A  B"1  =  A'  B1-1 
n  n  n  n 


then 


n 


SnA' 


n 


Bn  =  SnB^  n  =  0,1,2, 


(61) 


(62) 
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Proof : 


Substituting  equation  (62)  into  (59) 


SnA'  =  Sn_1K  A'  ,  +  SKVSn_1A'  ,  +  SSn“2A  0 
n  n  n-1  n  n-1  n-2 


and 


(63) 


SnB'  =  Sn-1K  B'  .  +  SK?  Sn_iB'  0 
n  n  n-1  n  n-2 


:n-lD, 


Dividing  equation  (63)  by  Sn,  recurrence  relation  (60)  is 
obtained.  (It  is  important  to  notice  the  generality  of 
equations  (59)  and  (60).  By  making  all  the  K's  =  0,  the 
recurrence  relationships  for  the  nth  convergents  for  the 
Cauer  first  form  is  obtained,  similarly  by  making  all  the 
K'is  =  0,  the  recurrence  relationships  for  the  nth  covergents 
for  the  Cauer  second  form  are  obtained.) 


and 


From  (59) 

X  =  [A  -  SIC  A  ,  -  SA  ,]A_1  ,  ' 
n  n  n  n-1  n-2  n-1 


K  =  [3  -  SK*  B  ,  -  SB  „ ]3~ 1  . 

n  n  n  n-1  n-2  n-1 


Solving  and  simplifying  terms 


A  B  ,  -  BA  ,  =  -S ( A  , B  „  -  B  .A  0) 
n  n-1  n  n-1  n-1  n-2  n-1  n-2 


so  by  induction 


A  B  -  BA  =  (-S)n-1*I 
n  n-1  n  n-i 


where  I  is  the  Identity  Matrix. 


(64) 


(65) 


(66) 
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Hence,  the  difference  between  two  consecutive  convergents  is 
given  by 


A  B-1  -  A  tB'1.  =  ( -S  )n~1[B1,B  .  ]-1 , 

n  n  n-1  n-1  K  K-l 


(67) 


Similarly,  from  (60),  and  following  the  same  procedure 


A'B'  ,  -  B'  A’  ,  =  i(A’  ,H  0  -  B'  ,A’  0 )  . 

n  n-1  n  n-1  s  n-1  n-2  n-1  n-2 


(68) 


So  by  induction 


A'  B'  .  -  3’  A'  , 
n  n-1  n  n-1 


i  (-  i)n-X 

s  s 


(69) 


and  the  difference  between  two  consecutive  convergen:s  is 
given  by 


A'  B71  -  A'  ,Bt1..  =  -  *  (-  i)n-1[3'  B’  ]_1 
n  n  n-1  n-1  s  s  n  n-1 


(70) 


By  looking  at  the  recurrence  relations  (59)  and  (60),  is 

observed  that  when  S=0  or  S==°,  respectively,  all  B^  and  all 

B'  are  non-zero  provided  that  all  K  and  K'  are  non-zero, 
n  r  n  n 

This  is  equivalent  to  apply  the  final  value  theorem  to 
equation  (67)  for  a  unit  step  inputs i  and  applying  the 
initial  value  theorem  to  equation  (70)  for  an  impulse 
response . 

Hence , 


ds 


t[A  B"1  -  A  , B_1, ]c_n  =  lim  —  I  *  S(-S)n-1[B  B  ,]-1' 
j  nr.  n-1  n-1  S  =  0  s  n  n-1 


s-»-0 

0 


j  =  0,1,...,  n-2 


(71) 
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Similarly , 


[A’  Bt1 
n  n 


A'  Bt1  1 
n-1  n-lJS=~ 


lim  S  *  i  I  (-J)Cff  E?  ,  D_1 / 
s  s  n  n-1  f 

S->oo  / 

0  j=0,l,..,n-2  J 

(72) 


Equations  (71)  and  (72)  show  the  results  that  were  expected, 
namely:  The  (n-l)th  approximant  and  the  nth  approximant 

goes  to  zero  in  the  sense  of  minimizing  the  follows  semi¬ 
norms  of  the  difference 


F  ,  (s)  I  |  , 

n-1  1 ' n-1 


Fn-l(s)Mn-l 


: 


n-2 
Z  |  — 
j=0  ds 


n-2 

Z 


3  Fn-l(s)ls=0 


3  =  0  ds 


-r-  F'  .  ( s )  _  =0 

3  n-1  's:“ 


(73) 


where 


and 


F  ,(s)  =AB-1-A  _B-1, 

n-1  n  n  n-1  n-1 


F*  -,  (  s  )  =  A'  B7 1  -  A’  B”1  , 
n-1  n  n  n-1  n-1 


(74) 


Note:  A  semi-norm  is  a  norm  which  does  not  satisfy  the  norm 

axiom  "||F(s)  ||  =  0  implies  F(s)  =  0". 

Since  the  derivatives  in  (73)  correspond  to  the  coeffi¬ 
cients  in  the  Taylor  series  expansion  of  the  functions 

F  .(s)  and  F’  ,(s)  about  the  points  S  =  0  and  S=°°.  Then,  it 
n-i  n-i 

can  be  deduced  that  the  output  difference  to  a  nth  approximant 
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has  the  same  Taylor  series  as  the  original  function  for 
terms  up  to  Sn~^  and  (^)n""\  effectively  neglecting  the  part 
of  the  transfer  function  which  differentiates  the  input  >_  n 
times  in  the  steady  state  as  well  as  in  the  transient.  This 
gives  the  sense  in  which  the  approximation  by  the  Cauer  mixed 
form  works . 


E.  SIMPLIFYING  A  MATRIX  TRANSFER  FUNCTION 

If  the  following  nth  order  system  is  given. 


T(s) 


CA2,nSn_1  + 

A  S  ^ 

. - -A24b 

*  A23S2 

+  A22S 

+  A21 

tAl,n+lSn  + 

A  S  ^ 

•  • -a14* 

+  A13s2 

+  A12S 

+  A11 

(75) 


a  simplified  model  of  the  system  is  desired.  By  performing 
the  Generalized  Routh's  Matrix  Array,  the  matrix  quotients 
of  the  Cauer  mixed  form  in  equation  (3)  can  be  obtained. 

If  an  m  order  is  desired  only  the  first  m  pairs  of  K's 
and  K'is  should  be  kept  in  equation  (3)  and  the  remaining 
should  be  omitted.  After  the  inverse  procedure  has  been 
performed  on  the  truncated  continuted  fraction,  the  simpli¬ 
fied  model  is  obtained. 

For  example  the  general  transfer  function  obtained  by  the 
matrix  continued  fraction  expansion  for  the  Cauer  mixed  form 
is  given  in  equation  (76). 

T(s)  =  [K1+K’1S+[K2|+K'2  +  [K3  +  K'3S+[K4i+K,4  +  [  .  .  .  D-1]-1]-1]"1]-1 
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If  a  second  order  simplified  model  is  required,  only  the 
first  two  pairs  of  matrix  quotients  will  be  kept,  that  is 


K  ,  K^,  1<2 ,  and  ,  and  the  rest  should  be  discarded. 

After  the  inverse  procedure  has  been  performed  on  the 
truncated  matrix  continued  fraction  expansion,  the  transfer 
function  given  in  equation  (77)  is  obtained.  Equation  (77) 
is  the  simplified  model  in  Mixed  Matrix  Cauer  Form. 

T(s)  =  [K'2S+K2]  x  [K’1K'2S2+(I+K1K'2+K'1K2)S+K1K2]''1  (77) 


The  truncation  procedure  outlined  above,  applies  in  similar 
manner  for  the  Cauer' s  first  and  Cauer 's  second  forms. 

This  methodology  is  particularly  advantageous  when  state 
space  terminology  is  used.  For  the  Cauer  mixed  form,  the 
state  space  formulation  is  written  as  follows: 


KiK2  K1K4 


^1^6 


K1K2  (Kl+K,3)K4  (K1+K3)K6 
K1K2  (K1+K3)K4  (K1+K3+K5)K6 


••KiKn 

.  .  (K*  +K'  )K’ 
ion 

.  .  (K’  +K’ +KV  )K' 
1  3  5  n 


KIK'  (K,,+K’)K!1  (K',+K'  +K'  )K*  .  .  .  (K*  +K’+K'  +  .  .  .+K'  ,  )K' 


12 

I+1<?1K2  +  K1K2  K1K4  +  K1K4 


1/ 

n 


.KIK  +K,  K* 
In  In 


K1K2+K1K2 

K1K2+K1K2 


K'1K2  +  K1K'2 


I+(K![+K'3)K4+(K1+K3)K'4, 

OWV'WH 


(K’1+K’3)K4*<K1*K3)K'4 


<K'*K',)K  +(K.+K,)K' 
ion  ion 

(KL+Kl, +KV)K  + 
loon 

(K.  +K,  +  Kj.  )K 

.  1  3  5  n 

(K' +K’+.  .  .+K'  ,  )K  + 
l  o  n-l  n 

( K,  +K,  +  .  .  .  +K  ,)K* 

1  3  n-l  n 


n/2 

E, 


Jn/2 


U2 


kxk2  kxk4 

K1K2  (K1+K3)K6 


.  . . K, K 
1  n 

. . . (K. +K, )K 
ion 


K1K2  <KltK3)K4  (WVK6  •■•<KltK3tK5)Kn 


KxK2  (K1*K3)K4  (KltK3tK5)K6  ...<K1*K3+...-MCn_1>Kn  E^J  |I 


Y  =  [K2,  K4,  K6, 


K  ]  E, 
n  3 


+  [K'2,  k^,  k;.. 


.  .  K’  ]  E0 
n  3 


where 


^n+ll 


•(n/2-l)m+l 

!(n/2-l)m+2 


'(n/2  )m 


I  =  identity  matrix,  and  E.,  ,  E^  ...  E^j*  E  ^ ^  2 ’  “  '  ^n/ 
have  the  same  dimensions  as  E^,  Ej,  •••  En/2  respectively 
and  m  corresponds  to  number  of  inputs  or  outputs. 


It  is  interesting  to  note  that  the  simplification  of  the 
state  equations  can  be  carried  out  by  partitioning  the  matrix, 
or  by  only  keeping  a  part  of  the  original  matrix  as  a  simpli¬ 


fied  model.  In  other  words,  if  a  simplified  model  of  a  two 
dimensional  matrix  is  required  which  is  equivalent  to  a 
2  x  (2m)  order  transfer  function,  the  upper  left  hand  corner 
of  the  original  matrix  is  taken  as  the  simplified  model. 


Therefore : 


KiK2 

KiK4 

5^2 

(K'+K'3)K 

K,K„ 

K.,K, 

1  2 

1  4 

K1K2 

(k1+k3)k 

I+K!lK2+K1K,2 

K1K2+K1K2 


I 

I 


U 


KiVKlK4 

I  +  (K'1+K,3)K4+(K1+K3)K,4 


(80) 


Y  =  CK2,  K4] 

where  the  vector  E  ,  E2,  E^,  E^,  E2,  keeping  the  same 

dimensions  as  were  given  for  the  equations  (78)  and  (79). 

The  method  used  for  the  simplifications  of  a  transfer  func¬ 
tion  and  state  space  formulation  based  on  the  Cauer  mixed 
form  holds  in  the  same  fashion  for  the  Cauer  first  and  second 
forms . 

For  effects  of  simulation  Figure  7  shows  a  diagram  for 
two  inputs  two  outputs  fourth  order  system  for  the  Cauer 
Matrix  mixed  form,  in  similar  form  for  the  Cauer  first  and 


+  Cic2 ,  K’4] 


(81) 
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Figure  7.  Simulation  diagram  for  fourth  order  system  with  m 
for  Cauer  matrix  mixed  form. 
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second  form  a  simulation  diagram  can  be  obtained  as  given 
in  Figures  8  and  9. 

It  is  important  to  notice  the  number  of  reductions,  q, 
which  one  can  realize  are  constrained  by  the  number  of  inputs, 
m,  and  the  order  of  the  system,  n,  where  q  =  n-km  for 
k  =  1,2,...  and  q  >  0 . 

E.  NUMERICAL  EXAMPLES  AND  GRAPHICAL  COMPARISON  OF  THE 
THREE  CAUER  FORMS 

1 .  Example  2.1 

Consider  two  inputs,  two  outputs  fourth  order  system 
where  a  reduced  second  order  system  is  required. 

The  state  space  equations  for  the  system  are: 


0 

0 

3.25 

L.125 

0 

0 

0 

0 

2.25 

-0.125 

0 

0 

X  = 

-1 

0.5 

-1.5 

-( 

) .  5 

X  + 

0.5 

0 

-1.5 

-0.5 

-1.0 

+  .0 

_ 

1 

0 . 5 

1 

-1 

1 

- 

1 

Y  = 

0 

1 

-1 

1 

X 

The  transfer  function  for  this  system  is  given  by  equation  (84). 


T(s)  = 

s  +  1 

s  - 1 

2s2  +  3s  +  2 

s  — 1 

2 

-2 

s+1 

-4s  -  4S  -  1 

2s  +  6s  +  3 

(84) 
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From  the  transfer  function,  the  matrix  Routh's  Array  and 
the  following  matrix  quotients  for  the  three  Cauer  forms  are 
obtained . 

Therefore,  reduced  second  order  models  are: 

Cauer  first  form 


1 

13 


Cauer  second  form 


2  1 

0.5  0 

= 

= 

_1  4_ 

2 

_0 . 12  5  0 . 2  5_ 

Cauer  mixed  form 


where  their  respective  transfer  functions  are  given  in  equa¬ 
tions  (84),  (85),  and  (86). 

Cauer  first  form: 


Fcl(s)  = 


l""3c  +  2 . 5  s  -  0.5 
s+0.5  s+2.5 


Cauer  second  form: 


2s  +  10s  +  6 . 5 


(84) 


Fc2(s) 


4s  +  2 
s  +  1 


-1 

2s  +  2 


9s  +  13s  +  5 


(85) 


» 
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Cauer  mixed  form 


F=3(s) 


3s  +  2 
s  +  1 


s  -  1 
s  +  2 


1 _ 

8s  +  5 


(36) 


Figure  10  presents  the  step  responses  of  reduced  models  just 
developed,  where  (10a)  corresponds  to  the  output  (1)  and  (10b) 
corresponds  to  the  output  (2). 

2 .  Example  2 . 2 

The  following  example  shows  the  reduced  systems  in 
Cauer  first  and  second  forms  are  only  stable  for  a  second 
order  approximation,  their  fourth  order  approximation  pro¬ 
vides  an  unstable  response.  Whereas,  the  fourth  order 
approximation  for  the  Cauer  mixed  form  is  completely  stable. 
Furthermore,  the  Cauer  first  form  gives  a  poor  approximation 
in  the  steady  state  portion  as  is  expected.  Consider  the 
following  sixth  order  system  with  two  inputs  two  outputs  as 
given  in  equations  (87)  and  (88). 


ro 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

1.0 

0 

0 

0 

0 

X  = 

0 

0 

0 

0 

1.0 

0 

x  + 

0 

0 

0 

0  • 

0 

0 

0 

1.0 

0 

0 

-50 

9.96 

-108.5 

19.6 

-16.5 

4 

2 

-1 

_4  9 

-19.96 

107 

-39 . 6 

13 

-10 

-2 

2 

10 

0.4 

1  0.2 

O 

O 

Y  = 

x . 

1.0 

4 

2  1 

1  0  _ 

(88) 
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Figure  10.  Outputs  for  the  original  fourth  order  system  and  reduced  second  order  system  in  the  three 
Cauer  forms. 


The  transfer  function  of  this  system  is  given  by  equation  (89). 


T  (  s )  = 


s  +  10 


0.2s  +  0.2s  +  0.1 


S  +  2s  +  1  s  +  4 


s3  +  10s2  +  55s  +  25  0.5s2  +  s2+0.2s  +0.2 

s3  +  3.5s2  +  1.5s  +  0.5  s3  +  6s2  +20  s  +  10 


-1 


or  if  the  inverse  is  to  performed 


(89) 


T(s) 


“0.25s5  +  0.1s4  +  14.9s3  +  79.25s2  +  209.75s  +  99.95 
_s5  +  7s4  +  25.5s3  +  40.5s2  +  33.5s  +  8 
0.2s5  +  1.7s4  +  7.1s3  +  6.8s2  +  8.48s  +  2.3 
-0.5s5  -  s4  +11. 3s3  +  93,58s2  +  244.76s  +  99.98 


0.5s6  +  13.25s3  +  130.55s4  +  562.53s3  +  1349 . 13s2+1049 . 87s+249 . 99 


(90) 


Cauer  Matrix  First  Form 


G 

G 


Performing  the  respective  Routh's  Array: 

;2] 


0 . 5 
1 


"  rio 

|_3.5 


1 

6 


[“55 

_  (_1.5 


0 

20 


0.2 

0 


]B  ][}• 

rs.5  -o.5  ”[  pa  -4.05 ~[  rss 

j_-3 . 5  4  J[_-49.5  15  ’  5  Jl°-! 


|“2 .686  -0.5  “j  |~9 .1 

[_-2.5  1.417j[_-3. 

E 


'7.823  -4.99S105 

-53.664  3 .153  J  (_?  .  5 


-3.444 
9  5  3.7  5 

— i 

0.02  “j 
10 


0 

4 

0 

10 


2] 

02~| 

J 


5  2 


* 


[7.2  34  -0.0  5  5") 
j_3.046  3.696J 

“25  0.02  "] 

_0 . 5  10  J 


and  from  the  Routh's  array,  the  respective  H’i's  are: 


2.5 

r 

0.035 

0 .054 

HI  = 

H'  = 

1 

_5 

2 

_0 . 1 9  7 

0 . 0  2  5_ 

2.566 

0.556 

0 .076 

-0.039 

H4  = 

_1 . 9  8  6 

3 . 526_ 

_-0  .  23 

0.006 

1.64 

-1.327" 

0 .29 

-0.006 

H6  = 

_-7 .73 

0.74  J 

_0.114 

0  .  369  _ 

The  state  space  representation  in  Cauer  first  form 
is  given  by  equations  (91)  and  (92). 


0 .284 

0  . 

16 

}  -0.09  -0.092 

1 

1 

0.338 

0.352") 

0.36 

0  . 

296 

i  -0.099  -0.189 

1 

1.562 

0.339 

0 .284 

'o'. 

16  " 

-0.052  -0.188 

1 

1.644 

0 . 544 

0.37 

0  . 

296 

-0.738  -0.246 

1 

1 

2.54 

1.629 

07235“' 

0. 

16 

-"oToS?-"-' 

57I3S“‘ 

1.967 

0 . 043 

_0. 37 

0  . 

296 

-0.738  - 

246 

0.387 

1.949 

p.  0  ! 

0 

0 

0  0 

ri  o~] 

0  1  ! 

0 

o  ! 

0  0 

0  1 

0  0 

1 

0  ! 

0  0 

z  + 

1  0 

U 

0  0 

0 

1  1 

0  0 

0  1 

0  0 

0 

0 

1  0 

I  5 

LP  0 

0 

0 

3  h 

LP  h 

(91) 
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y 


0.035  0.054 

0.197  0.024 


0,076  -0.039  |  0.289  -0.006 

I  z 

-0.28  0.006  |  0.114  0 . 369_ 

(92) 


The  reduced  fourth  and  second  order  model  are  obtained  by 
partition  of  the  respective  matrices  in  equations  (91)  and 
(92).  For  this  case,  the  fourth  order  model  is  unstable  and 
the  second  order  model  gives  an  untolerable  error  in  the 
steady  state  portion  of  the  response,  so  this  method  is  not 
applicable  for  this  example. 

Figure  11  shows  the  step  response  for  the  original 
system  and  the  reduced  second  order  model, 
b.  Cauer  Matrix  Second  Form 

Performing  the  respective  Routh ' s  Array: 


F25 
j_0.  5 

G# 

F52.61 
]_—  3.31 

F-0 .925 
[1.437 

F2 .663 
J_2  5.5 

F-0.175 
j_0 . 7  5  6 

G 


0.02  “ 
10 

F5  5 
|_1  •  5 

0.2  ~]  n.o 

20  J  L3  • 5 
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6 

0.1  " 
4 

G 
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0.7 

0 
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0 . 499~j  ri 
6.04  J  1__1 

0.7 

17.535 

1_0  .  9  9  5 

1 
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r-0.199 

0.096  1 

-0.397 

[_  0.738 

-0.245J 

5.1  " 
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0.5  ~j 

0.39 

0.127” 

-0.269 
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1  J 
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Figure  11.  Outputs  for  the  original  sixth  order  system  and  reduced  second  order  system  by  Cauer  first  form. 


and  from  the  Routh’s  array  the  respective  Hi's  are: 


2 . 505 

-0.058 

0.191 

0.008 

H1  = 

_-°  .  2 

2  .  5  0  5  _ 

H2  = 

_0 . 0  3  3 

0 . 229_ 

*-71 .026 

-9.128  “ 

“0.013 

-0 . 0  3 S' 

H3  = 

_-82 . 527 

-55.455^ 

H4  = 

_0 . 083 

0 . 065_ 

" 92.956 

2  5.02 

“0.6  04 

0.429" 

H5  = 

145.353 

67.357 

H6  = 

2.05 

-1.293 

The  state  space  representation  in  Cauer  second  form  is 
given  by  equations  (93)  and  (94). 


P0 .476 

-0.003 

|  -0.039 

0.098  i 

1.632 

-1.15 

[”1 

0] 

-0.045 

-0.571 

i  0.21 

-0.17  ! 

-5.255 

3.33 

0 

1 

r,  _ 

-0.476 

-0.008 

0.168 

-1.986  ; 

-22 . 573 

17.54 

l 

0 

-0.045 

-0.571 

-3.263 

0.322  ' 

58.535 

-32.97 

0 

1 

-0.476 

-0.008 

0.168 

-1.986 

-17.69 

10 

l 

0 

-0.045 

-0.571 

-3.263 

0.322 

3.3 

-8 .25 

0 

1 

(93) 

0.191  0.008  |  0.013  -0.038  j  -0.6  0.429-] 


0.033  0.229  |  -0.083  0.065  |  2.05  -1.293J 

(94) 


The  reduced  fourth  and  second  order  model  are 
obtained  by  partition  of  the  respective  matrices  in  equations 
(93)  and  (94).  For  this  case  like  the  Cauer  first  form,  the 
fourth  order  model  is  unstable  and  the  second  order  model 
offers  a  good  approximation  to  the  step  response  of  the 
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original  system.  Figure  12  shows  the  outputs  of  the  original 
and  reduced  second  order  model  for  a  step  input, 
c .  Cauer  Matrix  Mixed  Form 

Performing  the  respective  Routh ' s  Array: 


13-5  3  3 

•2]  G  s-2] 


[25  0.0  fj  [5  5  0-  21  flo  ri  n.  0.51 

[_0 . 5  10  J  LI-  5  20  J  |_3-5  6J  jjL  1  J 

G°  r]  B 

f  26.62  -4.49: 

L-54.31  13. o: 

["-1.832  —0.9  I 

L-7.2  -0.81J 

and  from  the  Routh' s  array  the  respective  Ki's  and  K’i's  are: 
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2 
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The  transfer  function  for  the  reduced  fourth  order  system 
is  given  in  equation  (95).  For  this  model  the  step  response 
is  completely  stable  and  the  approximation  has  an  error  of 
only  0.5%  over  the  entire  response  of  the  system.  Figure  13 


57 


Figure  12.  Outputs  for  the  original  sixth  and  reduced  second  order  system  by  Cauer  second  form. 


Figure  13.  Outputs  for  the  original  sixth  and  reduced  fourth  order  system  by  Cauer  mixed  form 


shows  the  outputs  of  the  original  and  reduced  fourth  order 
model  for  a  step  input. 


F(s) 


-0.012s3  -  0.063s2  +  1.44s  +  0.972 
_  0 . 061s3  +  0.773s2  +  0.278s  +  0.078 
0.012s3  +  0.173s2  +  0.068s  +  0.022" 
-0.03s3  -  0.234s2  +  1.784s  +  0.972 


0.03s 


0.9795s3  +  8.215s2  + 


S . 718s  +  2.431 


(95) 
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III.  DERIVATION’  OF  LOWER  ORDER  SYSTEM  FOR  THE  LINEAR 
REGULATOR  PROBLEM  VIA  CAUER  FORM 

A.  LINEAR  REGULATOR  SYSTEM 

R.  E.  Kalman  [12]  has  shown  that  when  a  system  is 
described  by  equations  (  96  )  with  a  performance  function 
given  by  equation  (  97  ) 

x  =  Ax  +  Bu  (  96  ) 

J  =  15l[xTQx  +  uTRu]dt  (97) 

and  if  the  plant  is  completely  controllable  then  an  optimal 
u*  exists  and  is  given  by  equation  (  98  ) 

u*  =  -Kx  (  98  ) 

where  the  feedback  matrix  X  for  the  controller  depicted  in 
Figure  14  is  a  constant  matrix  as  t 

In  recent  papers,  Goldman  [10]  and  Aoki  [18]  have 
derived  a  way  to  obtain  near  optimal  solutions  using  the 
"aggregation  matrix" ,  to  transfer  an  arbitrary  state  x  to  the 
origin  of  the  state  space  while  minimizing  the  criterion 
function  given  in  equation  (  97  ) . 

The  present  work  will  show  that  the  reduced  order  optimal 
regulator  can  be  obtained  where  the  original  system  is 
translated  in  Cauer  form  and  also  a  near  optimal  solution 
can  be  found  through  a  mere  partition  of  its  optimal  solution. 
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Figure  14  .  Linear  Regulator  System 

3.  THE  STATE  VARIABLE  REPRESENTATION  BY  CAUER 
SECOND  FORM  AND  ITS  OPTIMAL  SOLUTION 

It  is  known  that  the  standard  state  variable  representa¬ 
tion  of  a  linear  time-invariant  system  is  given  by  equations 
(  99  )  and  (100  )  . 

x  =  Ax  +  Bu  (  99  ) 

y  =  Cx  (100) 

T 

where  x  is  a  column  vector  nxl  given  by  Xj.-.x^]  ,  A 

is  an  nxn  constant  coefficients  matrix,  x  is  a  column  vector 
nxl  given  by  [x^,  X2>--*x  3  ,  B  is  an  nxl  constant  coefficients 


Goldman  (10),  Chin  and  Shieh  (11)  have  proven  that  the 
system  given  by  equations  (  99  )  and  (100)  can  be  represented 
in  Cauer  second  form  as  shown  in  equation  (101)  and  (102). 


2  = 

Fz  +  Gu 

(101  ) 

y  = 

Mz 

(102  ) 

where , 

F  = 

PAP"1 

(103  ) 

G  = 

PB 

(104  ) 

M  = 

CP"1 

(135  ) 

where  the  matrix  P  is  an  nxn  upper  triangular  matrix  and  the 
elements  in  the  triangle  are  copied  directly  from  the  elements 
of  the  Routh's  array,  where  the  nth  row  of  the  P  matrix 
is  the  (2n+l)th  row  of  the  Routh's  array  and  the  states 
variables  x  and  z  are  related  by  equation  (106). 


V  ‘ 


f 


V- 

i  I 


2  =  Px 


(106) 


If  the  system  in  ( 99  )  and  (100)  is  put  in  transfer  function 
notation  the  result  is  given  by  ( l Q7  >  - 


_a21+a22S+a23S2+a24S3+.  .  .+a2nSn~1 
all+al2S+al3S2+al4S3+. . .+a2n+lSn 


(107) 


■» 


_/ 


J 
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The  Routh 1 s  array  can  be  formed  as 


all 

al2 

al3 

. . .  aln 

a21 

a22 

a23 

.  .  .  a2n 

(  108) 


and  the  elements  of  third,  fourth,  and  subsequent  rows  can 
be  evaluated  from  the  following  algorithm. 


ajk  =  aj-2,  k+l-f j -2aj -1 ,  k+1  j = 3 ,4 , . . . ,n+l 

k=l ,  2  ,  .  .  . 


fp 


ap,l 
ap+1 ,1 


p= 1 , 2  , .  .  .  n 


(10  9  ) 

(110  ) 


The  elements  of  the  matrix  F  can  be  obtained  also  from 
equation  (111)  where  ap,l  and  ap+1,1  correspond  to  the 
elements  of  the  first  column  of  the  Routh's  array  and  then 
F  is  formed  as  shown  in  equation  (111). 


F=- 


”f  2fl 
f  2f  1 
f  2fl 


f  4f  1 

f 4 ( f  1+f  3 ) 
f 4 ( f 1  +  f  3 ) 


f 6f 1  ...fnfl 

f 6(fl+f 3)  . . . fn( f 1+f 3 ) 

f 6 ( f 1+f  3  +  f 5 ) . . . fn( f 1  +  f 3  +  f 5 ) 


f 2f 1  f 4 ( f 1+f 3 )  f6(fl+f3+f5) . . . f n( f 1+f 3+ . . . fn-1 ) 


(111) 


The  elements  of  the  matrix  G  correspond  to  a  column  vector 
of  nxl  of  l's  and  the  elements  of  the  matrix  M  can  be  formed 
also  by  fp’s  where  p=2 ,4 ,6 , . . . 2n  as  shown  in  equation  (112). 


M  =  [f 2 ,  f 4 ,  f 6 , . . . f 2n] 


(112  ) 
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The  performance  function  given  in  equation  (  97  )  can  be 
translated  to  Cauer  form  by  using  the  relationship  given 
in  equation  (106),  i.e., 


J  =  Q  z*  +  u*T  RU*  ]dt 

c  c 

where 

Q  =  Q 

and  the  optimal  u*  is  given  by  equation  (115  )  . 

u*  =  -K  z* 
c 

where 

K  =  K  P-1 
c 

So  the  equivalent  system  in  Cauer  form  given  in  Figure  14 
can  be  represented  by  Figure  15 . 

|  "plant 


| REGULATOR  IN_  CAUER _ 1_ 


Figure  15 .  Linear  Regulator  in  Cauer  Form 


Note  that  Kc  is  the  feedback  gain  matrix  of  the  Cauer 
system  and  u*  is  the  optimal  control  law  of  both  systems 
(Figure  14  and  15). 

From  the  theory  of  optimal  control,  z*  and  u*  are  found 
by  solving  the  set  of  necessary  conditions: 

Z*  =  ~r~  (z* ,  u* ,  p* ,  t) 

3p* 


0 


0 


3  H 


( z*,  u*,  p*,  t) 

(z*,  u*,  p*,  t) 


for  all 
te[fO ,ff 3 


(117  ) 


Where  the  function  H  is  called  the  Hamiltonion  and  is 
defined  by: 

H  (  z  ,u  ,p  ,  t )  =JsCz^Qcz+U^RU]  +  p^[Fz  +  Gu] 

(118) 

and  p  is  the  Lagrange  multiplier  or  Costate  state. 

Substituting  equation  (118)  into  equation  (117),  the 
canonical  system  in  equation  (119)  is  formed. 


z* 

F  -GR_1GT 

z* 

p* 

-Q  -FT 

^  c  -1 

-  p*  _ 

(119  ) 
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Fortunately,  for  the  Optimal  Regulator  problem,  it  is  not 
necessary  to  solve  these  equations.  Kalman  [12],  Tyler, 
J.S.  and  Tuteur,  F.B.  [19]  have  shown  that  when  the  optimal 
control  u*  is  generated  by  equation  (115),  the  solution 
for  Kc  is  obtained  from  equation  (120)  where  the  matrix,  T, 
is  the  solution  to  equation  (121)  in  steady  state. 


K  =  R_1GTT  (120) 

G 

0  =  ftt+tf-tgr_1gtt+q  (121) 

c 


Equation  (121 )  is  the  steady  state  form  of  the  Ricatti 
equation . 

For  the  single  input  case  the  weighting  matrix,  R,  is 
a  scalar.  When  this  is  the  case,  a  lower  order  optimal 
linear  regulator  can  be  found  just  by  mere  partition  of 
F,  G,  M,  Qc  and  Kc  matrices  as  shown  in  equation  (122). 
Hence  the  new  system  is  depicted  in  Figure  16  . 


Figure  15  .  Reduced  Linear  Regulator  in  Cauer  Form 
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M  =  [ f 2  f 4  :  f 6 . . .f2n] 

qcll 

qcl2 

!  ‘ • ' qcln 

I 

V 

qc21 

qc22 

!  qc2n 

i 

K=[k,,  k  ;  k  ,  ] 

c  ell  cl2  ,  cln 

qcnl 

qcn2 

•  •  a  Q 

xnn 

(122  ) 

The  reduction  scheme  outline  above  has  left  u  *  and  u  * 

c  r 

unchanged.  Therefore;  u  *  is  the  suboptimal  solution  for 
all  three  systems,  that  is  the  original  system,  the  trans¬ 
formed  system  and  the  reduced  system. 

C.  THE  REDUCED  ORDER  SOLUTION 

Based  on  equation  (122),  the  reduced  system  is  given  by 

zr  =  Frz  +  Gr  u  (123) 


y  =  Mrz 


(124) 


where  Fr,  Gr  and  Mr  are  the  partitions  of  the  matrices  F,  G 
and  M  respectively. 


68 


Goldman  [10]  has  shown  that  the  original  system  is 
related  to  the  reduced  order  model  by  the  following  equations: 

zr  =  eAe+zr  +  eBu  (125  ) 


y  =  Ce  zr 


x  =  e  zr 


zr  =  ex 


(126  ) 
(127  ) 

(12  8  ) 


where  the  matrices  e  and  e+  are  partition  of  the  P  and  P  ^ 
matrices  in  rectangular  form  respectively,  such  that  e  is  of 
order  (rxn)  and  e+  is  (nxr) ,  where,  n,  is  the  order  of  the 
original  system  and,  r,  is  the  order  of  the  reduced  system. 

From  equation  (127  )  and  (  98  )  it  is  easy  to  show  that 

u*  =  Ke*  zr  (12  9  ) 

then 

Kr  =  Ke+  (130  ) 

or 

K  =  Xre  (131  ) 


Note  that  from  equations  (98  )  and  (106  ) 

K  =  KP'1  (132) 

c 

or 

K  =  K  P  (133) 

c 
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Then  from  equations  (130)  and  (133) 


Kr  =  Kc  Ir  (134) 

where 

Ir  =  Pe+  (135) 

The  matrix  Ir  has  dimension  (nxr)  and  is  a  special  type  of 
identity  matrix  where  the  last  (n-r)  rows  are  zeros  and  the 
first  rxr  components  correspond  to  the  identity  matrix.  In 
the  same  way,  it  is  possible  to  show  that  the  different 
relationships  that  exist  for  the  weighting  matrices,  Q,  Qc 
and  Qr  are  given  in  equations  (136),  (137),  (138)  and  (139). 


Qc  = 

P"1T  Q  P_1 

(136) 

Qr  = 

e+TQe+ 

(137) 

Qr  = 

+T„T„  „  + 
e  P  QcPe 

(138) 

Qr  = 

Ir^  Qc  Ir 

(139) 

1 .  Example  3 . 1 

The  transfer  function  of  a  system  is  given  in  equation 
(l40 ) .  It  is  desired  to  control  the  system  in  such  a  way  that 
a  performance  function  given  in  equation  (141)  is  minimized. 

Due  to  the  systems  complexity,  a  near  optimal  solution  is 
desired . 

H(s)  =  \ -  (140) 

SJ  +  6  .IS  +5. 6S  +  0 . 5 
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(141) 


J  = 

/1sCxTQx+UTRU]dt 

where 

"5 

0 

0" 

Q  = 

0 

4 

0 

R  = 

1 

_0 

0 

1_ 

From  the  transfer  function 

ro 

1 

0 

"  0  “ 

A* 

0 

0 

1 

B= 

0 

L-o 

.  5 

- 

5.6  -6.1. 

_  1  _ 

The  Routh 

'  s 

array 

0. 

5 

5 

.  6 

i — 1 

10 

1.0 

1.0 


0 .0 


{5 . 6 


6.1 


0.0 

TTol 


-1.09  -0.179 

gTl82  l.'0| 


0.932 

fTTo] 


From  the  Routh 

s  array 

5 . 6 

6.1 

l.(T 

P  = 

0 

5.182 

1.0 

0 

0 

i.o_ 

C=  [1  0  0] 


(142  > 


I 


p 


-1 


0.179 

0 

0 


-0 . 2102 

0.193 

0 


0.032 

-0.192 

1.0 


then 


e  - 


5.6 

0 


6.1 

5.182 


1.0 

1.0 


fo  .17  9 


0 


-0.2102 
0.19  3 
0 


The  Qc  matrix  obtained  from  equation  (136)  is 


"0.159 

-0.188  ; 

0.029  ” 

Qc  = 

-0.188 

0.372  1 

.  i 

-0.183 

0.029 

-0.183 

1.154 

The  Fc,  Gc  and  Me  matrices  from  equations  (103),  (104)  and 
( 10  5 )  are 


10.089 

0.105 

i 

i 

0.015  “ 

Fc  = 

-0.089 

-0.978 

1 

1 

0.139 

_-0 .089 

-0.978 

-5 . 031_ 

Gc  = 

~  1  “ 

1 

:r_ 

Mc=  [0. 

1786 

-0.211 
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(143  ) 


(144) 


(145) 


(146) 


(147) 


» 


p 


w 


The  reduced  system  is  given  by  equation  (148)  and  the 
corresponding  matrices  that  describe  the  system  in  equation 
(149) . 

TTll/  -0 . 0316s+0 . 19  3 

H4(s)  =  - - 

s  +1.065S+0.0965  (148) 


-0.089 


Fr  = 


-0.089 


0.105“ 

-0.978 


1 


Gr  = 


1 


Mr=  [0.1786-0.211] 


0.159 


Qr= 


-0.188 


-0.188 

0.372 


(149  ) 


The  computer  solution  for  the  three  systems  for  the  feedback 
gain  are  shown  in  Table  I. 


TABLE  I 


SYSTEMS 

OPTIMAL  FEEDBACK  MATRIX  K 

ORIGINAL 

1.79128425 

2.087266 

0.41034 

CAUER  FORM 

0 . 319649369 

0.0243818 

0.06619982 

REDUCED  CAUER  FORM 

0.31965241 

0.02258695 
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Table  II  shows  values  found  from  the  reduced  system 
back  to  Cauer  form  and  original  system  with  the  equation 


(131)  and  (132). 


TABLE  II 


SYSTEMS 

NEAR  OPTIMAL  FEEDBACK  MATRIX 

ORIGINAL 

1.7908 

2.067 

0.342 

CAUER  FORM 

0.31965 

0.022586 

0.0 

REDUCED  CAUER  FORM 

0.31965 

0.022586 

Note  that  a  near  optimal  feedback  matrix  is  obtained  directly 
from  equation  (131)  but  not  to  the  Cauer  system  since  the 
reduction  technique  discarded  the  last  terms . 

Figure  17  shows  the  step  response  for  the  optimal  and 
reduced  linear  regulator.  Figure  18  shows  the  optimal  and 
suboptimal  control  laws .  Figure  19  shows  the  step  response 
for  the  optimal  and  suboptimal  linear  regulator. 

2 .  Example  3 . 2 

The  methodology  just  developed  is  applied  to  a  simple 
model  of  a  nuclear  reactor  power  generator  [20,  21]. 

The  heat  generating  process  of  a  nuclear  reactor  is 
dependent  upon  the  mechanism  called  fission  (a  fragmentation 
of  matter).  The  power  generated  by  this  process  is  directly 
related  to  the  population  of  neutrons,  n(t)  and  can  be 
described  by  the  following  differential  equation  (developed 
from  a  diffusion  balance  equation) . 


)  n(t)  +  Xc(t) 


(150  ) 


n(t)  = 

c ( t )  =  ( j)n(t)  -  Xc( t )  (151) 

where 

fik(t)  =  <5 kc(t)-an(t)  (152) 

The  variable  <5k(t)  is  the  input  to  the  process  and  is  given 
the  name  "reactivity1'.  It  is  clear  by  inspection  that 
5k(t)<6  for  stable  system  (in  a  linear  sense).  The 
variable  c(t)  is  a  measure  of  the  concentration  of  fragments 
(precursors)  that  produce  delayed  neutrons  according  to  a 
time  delay  (1/X)  called  the  "half-life"  of  the  precursor. 

The  input,  5kc(t),  is  the  control  (reactivity)  that  is 
associated  with  the  control  rod  position  and,  a,  is  a 
temperature  feedback  (reactivity)  coefficient. 

The  parameters  for  this  system  are: 
a  =  IQ"5  kw-1 
1  =  10~3  sec. 

8  =  0.0065 
X  =  0.1  sec."3- 

at  t=0  n(0)=10kw  (as  the  operating  output  in  steady 

state  conditions) 

The  control  problem  is  stated  as: 

Find  the  Optimal  Control  Policy  u*(t)  that  will 
transfer  the  power  level  n(t)  from  the  operating  level 
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a 


n(0)  =  lOkw  to  a  new  level  n(t)  =  50  kw  where 


u(t )  =  5  kc(t) 


(153) 


using  the  performance  measure  given  in  equation  (  97  )  where 


(154) 


R  =  1 


and  compare  this  solution  to  the  near  optimal  one  obtained 
from  the  reduced  order  model. 

From  the  state  of  the  problem  the  following  diagram 


is  drawn 


Figure  20.  Block  Diagram  Optimal  Control  of  a  Nuclear 
Reactor . 
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The  double  subscript,  ss,  means  steady  state  and  the 
new  definition  of  state  variables  are: 

xl(t)  =  n(t)  -  n  (155) 

ss 

x2(t)  -  c(t)  -  c__  (156) 

ss 

x3(t)  =  5  kc(t)  -  5 kc  '  (157) 

s  s 

The  steady  state  conditions  are 

n(t )  =  0  (158) 

c ( t )  =  0  (159) 

5  k( t )  =  0  (160) 

The  initial  conditions  at  n(t)  =  50kw  are: 

from  the  state  variables  definitions,  equations  (155), 

(156)  and  (157)  and  equations  (150),  (151)  and  (152). 

c(t)  =  ~  n(t)  -  c(t)  (161) 

c( 0 )  =  n( 0 )  (162) 

c( 0 )  =  640  (163) 

x2 ( 0 )  =  c(0)  -  c  (164) 

S  5 

x2 ( 0 )  =  -2560  (165) 

i  k(t)  =  an(t)  -  A1  +  8  (166) 

6  k(  0 )  =  an(  0 )  -  XI  6  (167) 
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•*  «. ; 


6k(Q) 

=  10~4 

5  k 

=  dn 

ss 

ss 

5  k 

ss 

=  5x10 

ss 


ss 


+  3 


-4 


x3 ( 0 )  =  5  k(0)  -5k 


ss 


x3(0)  =  -4x10 

xlCO)  =  n(0)-n  _ 
ss 

xl( 0)  =  -40 


then 


xl(0f 

~-40 

x(0)  = 

x2(  0) 

= 

-2560 

x3(0)_ 

-u 

-  4x10 

and  the  optimal  control  law  U*(t)  will  be 


u*(t)  =  -Kx(t) 


(168) 

(169) 

(170) 

(171) 

(172) 

(173) 

(174) 


(175) 


same  as  given  in  equation  (  98  ).  The  system  in  variable 
form  is  given  by  equations  (  99  )  and  (100).  Since  equations 
(150)  through  (157)  are  non-linear  linearization  is  required 
(see  appendix  A).  The  linearized  system  is  represented  by 
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0 


-  kan  +8  ) 

1  ss 

A 

n 

ss 

1 

x(  t  )  = 

6  /I 

-A 

0 

x  + 

0 

0 

0 

(176  ) 


y(t)  =  Cl  0  0]x(t) 


where 


-6.9 


A 


6.4 


0 


0.1  5xl04 

-0.1  0 
0  0 


(  177) 


(  178) 


B  = 


(179) 


C  =  [1  0  0] 


(180) 


Solving  the  Ricatti  equations,  the  matrix  K  is  given  by- 
equation  (131). 

K  =  [2.5xl0-4  8xl0"6  7.05]  (181) 

The  transfer  function  for  the  optimal  system  will  be: 


H(s) 


7S 


5x10  (S+0.1) 

7T, 


S  +14.05S  +61.9S+4.16 


(182) 
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The  Routh's  array  for  this  system  is: 


4.16 

61.9 

14.05 

5000 

50000 

0.0 

pr: 

3 

14.05 

1.0| 

46539 

-246.3 

II4.T57 - 

1.0 

-3533 . 58 

Fo| 

From  the  Routh' 

'  s  array 

“20 . 3 

14.05 

1.0" 

P  = 

0.0 

14.157 

1.0 

_0.0 

0.0 

1.0. 

“0.04  9 

1 

-0.049 

-0.0004" 

P"1  = 

0.0 

0.07 

-0.07 

0.0 

0.0 

1.0 

(183) 


(184  ) 


and  the  rectangular  partition  of  P  and  P  ^  matrices  give 


20.3 

14.05 

1.0“ 

_0.0 

14.157 

1.0 

(185) 

”0.049 

-0.049" 

0.0 

0.07 

0.0 

0.0 

(186  ) 
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and  from  equations  (125)  and  (126)  the  reduced  second  order 
suboptimal  system  will  be 


-0.205 

-2.735" 

1 

zr(t)  = 

zr  ( t )  + 

_-0 .205 

-4 . 169_ 

_  1  _ 

y ( t)  =  [246.305  3287.5]zr(t) 

The  corresponding  transfer  function  is  given  in  equation 
(189  >  • 


„  ,  .  3533.58S+353.17 

Hrop ( 3 ) =— * - 

S  +4. 3738S+0 . 294 

In  Figure  21  is  shown  the  step  response  for  the  original 
optimal  and  reduced  suboptimal  system. 
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(187  ) 

(188  ) 

(189  ) 


.00 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


Comparative  analysis  of  Cauer  forms  methodology  for 
multivariable  system  reduction  is  established.  The  developed 
methodology  is  based  on  the  Cauer  matrix  generalized  form, 
which  offers  the  closest  approximation  to  the  original 
system.  The  proposed  methodology  and  because  of  the  nature 
of  the  Cauer  mixed  form  shows  it  to  be  superior  to  any  other 
method  proposed  to  date  since  it  provides  satisfactory 
results  for  both  the  transient  and  the  steady  state  portion 
of  the  system  response.  The  methodologies  in  state  space 
as  well  as  in  the  S  domain  are  developed.  For  the  basis 
of  comparison,  reduced  order  models  using  the  three  Cauers 
forms  are  obtained  for  two  different  examples.  The  results 
clearly  show  the  superiority  of  the  Cauer  mixed  form  over 
the  entire  frequency  range  of  system's  response.  The 
proposed  methodology  is  algorithmic  therefore,  it  is 
amenable  to  digital  computation. 

A  lower  order  optimal  linear  regulator  can  be  obtained 
by  mere  translation  of  the  original  system  to  Cauer  second 
form  and  their  partition  of  the  different  matrices  as  shown 
in  Section  III-B.  A  suboptimal  feedback  matrix  for  the 
original  system  can  be  obtained  by  multiplication  of  the 
lower  feedback  optimal  matrix  by  the  rectangular  partition 
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of  the  matrix  P.  The  responses  of  both  systems,  the  original 
and  the  reduced,  as  well  as  optimal  and  sub-optimal  systems 
are  extremely  close. 


The  results  presented  here  are  encouraging,  there  is 
a  distinct  need  for  future  research,  in  particular  in  the 
reduction  of  systems  with  any  given  number  of  inputs  or 
outputs . 


APPENDIX  A 


Both  functions,  F  and  g  can  be  expanded  in  a  Taylor  series 

about  the  points  (x  ,  u  )  resulting  in  the  following 

s  s  ss 

representation  of  the  system  equations. 
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x  +  <$x=F(x  ,  u  )+A5  x+gS  u+a(S x,<S  u) 

ss  ss  ss 

y  +<$y=g(x  ,  u  )+GS  x+D5  u+B  (6  x  ,6  u) 

3  S  S  S  S  S 

where : 

A  =  nxn  matrix 
B  =  nxm  matrix 
C  =  vxn  matrix 
D  =  vxm  matrix 

n  =  order  of  the  system,  m  =  number  on  inputs,  v=number  of 
outputs. 

The  functions  a(6x,du)  andg(Sx,6u)  represents  all 

second  order  and  higher  terms  in  the  Taylor  series  expansion. 

Substitution  of  equations  (A-l)  and  (A-2)  evaluated  at 

(x  ,  u  )  and  neglecting  second  order  terms  and  higher 
s  s  s  s 

terms,  yields  the  following  perturbed  equations  of  motion 


5  x  =  A5  x  + 

B6  u 

( A-9 ) 

6  y  =  OS  x  + 

D5  u 

(A-10) 

The  two  equations  above  approximate  the  dynamic  behavior 

of  the  system  about  the  operating  point  (x  ,  u  ) .  The 

s  s  s  s 

elements  of  the  motion  A,  B,  C  and  D  are  given  by 

Al3  =  fxj  1  =  1,2'*,n  S  =  1,2 ' * ,n 


(A- 7) 

(A-8  ) 
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